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ABSTRACT 

In this paper we analyse different chronometers based on the models of chemical 
evolution developed in Chamcham, Pitts & Tayler (1993; hereafter CPT) and Cham- 
cham & Tayler (1994; hereafter CT). In those papers we discussed the ability of our 
models to reproduce the observed G-dwarf distribution in the solar neighbourhood, 
age-metallicity relation and radial chemical abundance gradients. We now examine 
their response to the predictions of cosmochronology. We use the recent production 
ratios of the actinide pairs 235 TJ/ 238 U and 232 Th/ 238 U provided by Cowan, Thiele- 
mann & Truran (1991) and the observed abundance ratios from Anders & Grevesse 
(1989) to determine the duration of nucleosynthesis in the solar neighbourhood, and 
thus to determine maximum likelihood estimates and confidence intervals for the infall 
parameter, /3, which controls the growth rate of the disc in our models. We compare 
our predictions for the age of the disc with the age of the galaxy estimated from mod- 
els of white dwarf cooling and from the age of globular clusters. From our statistical 
analysis we find that these three methods of age prediction appear to be consistent for 
a range of maximum likelihood values of [3 which is in good agreement with the values 
considered in CPT and CT, which were found to give a good fit to the observational 
data examined in those papers. We also briefly consider the consistency of our results 
with the age of the universe predicted in different cosmological models — a topic which 
we will investigate more fully in future work. 
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1 INTRODUCTION 

In a series of earlier papers, models for the chemical evo- 
lution of the galaxy have been developed and studied by 
Chamcham, Pitts and Tayler based on the adoption of a 
threshhold - derived from the Toomre stability criterion - 
for the onset of star formation (SF) in the galactic disc. 
In Chamcham, Pitts & Tayler (1993; hereafter CPT) and 
Chamcham & Tayler (1994; hereafter CT) the ability of 
these models to reproduce the observed G-dwarf distribu- 
tion in the solar neighbourhood, age-metallicity relation and 
radial chemical abundance gradients was discussed. In this 
paper we examine the response of our models to the predic- 
tions of cosmochronology for the age of the galactic disc. 

The age of the galaxy continues to be one of the most 
debated subjects in astronomy - particularly in the light of a 
number of recent determinations of the Hubble constant (cf. 
Schmidt, Kirshner & Eastman 1994; Freedman et al. 1994) 
which appear to indicate an age of the universe of less than 
10 Gyr for the standard inflationary model of a flat universe 



with zero cosmological constant. These results are in stark 
contradiction with astrophysical constraints on the age of 
the galaxy from, e.g., the cooling of white dwarfs and the 
age of globular clusters (cf. Renzini 1993). If both age de- 
terminations were correct then one would reach the absurd 
conclusion that the universe is younger than our galaxy, and 
the resolution of this paradox is one of the outstanding prob- 
lems of observational cosmology today. 

In Table 1 we summarise the constraints on the age of 
the universe reported in recent literature, as provided by the 
four methods of age determination which we consider in this 
paper. Of course, of these four methods, only by determining 
the values of cosmological parameters can one directly infer 
the age of the universe - and even then only in the context of 
a given cosmological model. As indicated in the first column 
of Table 1, the other three methods constrain the age of the 
solar neighbourhood, galactic disc and galaxy respectively 
- each of which may then be regarded as a (progressively 
larger) lower bound for the age of the universe. Thus, al- 
though the four methods are not directly comparable, one 
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can at least demand that their predictions should be phys- 
ically consistent with each other - a useful approach which 
was adopted in, e.g., Tayler (1986) and which we also adopt 
in this paper. It would obviously be unreasonable for the 
age of the galaxy to be smaller than the age of the solar 
neighbourhood, or larger than the age of the universe, for 
example. 

With respect to the second and fourth rows of Table 
1, it would also be naive to imagine that all values in the 
cited range are equally likely. The age determined from cos- 
mological models extends to the limits of the quoted range 
only by a somewhat delicate 'balancing act' between the 
values of the relevant cosmological parameters. In particu- 
lar, if the Hubble constant is indeed high - say, Ho > 75 - 
then all but the extreme lower limit of the quoted range is 
strongly excluded in the standard inflationary model. Simi- 
larly, if one believes the measured values and observational 
errors of the actinide pairs' abundance ratios which we con- 
sider in this paper, then we shall see that a careful statistical 
analysis can place tighter constraints on the most likely age 
range determined from cosmochronology - at least within 
the framework of a given chemical evolution model - than 
suggested by the full extent of the range quoted in Table 1. 

In what follows we shall focus on the method of cos- 
mochronology and compare its results with the predictions 
of the other methods. Cosmochronology has recently been 
gaining in credibility, especially with its extension to indi- 
vidual stars by Butcher (1987) - although applications have 
usually been limited to the solar system - and the progress 
of theoretical and experimental nuclear physics (Klapdor- 
Kleingrothaus 1991). This method is, however, still strongly 
dependent upon one's choice of chemical evolution model. 

In CPT and CT we have shown that one of the weak- 
nesses of our models is the number of free parameters which 
we are using. The aim of this paper is, therefore, to place 
some constraints on the most important of these parameters 
through the analysis of cosmochronometers - specifically, 
the abundance ratios of two actinide pairs - by confronting 
their predicted age of the Galactic disc with other astrophys- 
ical and cosmological predictions. The particular strength of 
our models, on the other hand, is that the star formation 
(SF) law is deduced consistently with some physical input, 
instead of the simple Schmidt law (e.g. SF proportional to 
the gas density) commonly used by different authors in their 
cosmochronological models (Schramm & Wasserburg 1970; 
Clayton 1988; Malaney & Fowler 1987; Pagel 1989). 

Our present approach is analogous to that of Mathews 
and Schramm (1993) in that both models consider a delay in 
SF instead of an instantaneous birth of stars in the Galaxy. 
In our models this delay is related to the stability properties 
of the disc (cf. CPT), whereas in the model of Mathews 
& Schramm it is introduced only schematically to sketch 
the merger process in the formation of the Galaxy. In our 
models this delay varies between 1-5 Gyr at the solar 
neighbourhood, depending on the growth time of the disc 
(or equivalently the rate of infall of material onto the disc) 
and the relative roles of heating and cooling processes which 
regulate the star formation rate (SFR). 

The structure of this paper is as follows. In section 2 we 
describe briefly our present chemical evolution model, sum- 
marising the relevant details from earlier papers. In section 
3 we discuss the results of applying our chemical evolution 



model to predict the abundance ratios of the two chronome- 
ters as a function of the age of the disc, and compare these 
model predictions with the observed abundance ratios. We 
then describe a method for obtaining maximum likelihood 
estimates and confidence intervals for the infall parameter, 
/3, which controls the growth rate of the disc in our models, 
as a function of the age of the disc. In section 4 we review 
in more detail the constraints on the age of the galaxy pro- 
vided from white dwarfs and globular clusters. In section 5 
we then discuss the consistency of the age predictions from 
our models with those from white dwarfs and globular clus- 
ters, and on this basis consider what contraints are placed 
on the parameter /3, and on the age of the disc. We also 
briefly discuss the consistency of our results with the age of 
the universe predicted from the values of cosmological pa- 
rameters - a subject which we will consider in more detail 
in a separate paper. Finally in section 6 we summarise our 
conclusions. 



2 CHEMICAL EVOLUTION MODEL 
2.1 Threshold of star formation 

The determination of the age of the Galaxy by the method of 
cosmochronology provides at least a lower limit for the age of 
the universe. This method is based on the assumption that 
the earth and meteorites solidified some 4.5 Gyr ago and 
that the isolation and condensation of the solar system from 
the ISM occurred at Tq = 4.6±0.1 Gyr ago. The present-day 
abundance ratios of the chronometric pairs 235 u/ 238 U and 
232 Th/ 238 U are determined from the analysis of meteorites 
(i.e. carbonaceous chondrites) and moon rock samples. Their 
values at the time of the formation of the solar system are 
extrapolated back in time and are respectively 0.33 ± 0.03 
and 2.32 ±0.23 (Anders & Grevesse 1989). 

The duration of nucleosynthesis, t n uc, from the start of 
SF in the Galaxy until the formation of the solar system is 
determined by a suitable chemical evolution model. One can 
estimate from this a lower limit for the age of the galactic 
disc to be 

Td = Tq + t nuc (1) 

More correctly, equation [l] accounts for the age of the solar 
neighbourhood - in other words the time since SF began at 
our location in the Galaxy. We thus rename this estimated 
age T s instead of T D . Q 

The variation in the models of cosmochronology used by 
different authors stems from the different relations assumed 
between the SFR and the gas mass, and their choice of model 
for the infall of material onto the disc (Clayton 1988). To 
our knowledge all previous authors have assumed that SF 
begins either instantaneously everywhere in the Galaxy or 
in the case of Mathews & Schramm (1993) - after an 
arbitrary delay of a few Gyr. As shown in the threshold 
models of CPT, however, the delay, 8, in the onset of SF 
in the solar neighbourhood is not an arbitrary parameter 
but is determined self-consistently by the other parameters 
of our chemical evolution model. We can therefore deduce a 

* Note that Tg is therefore comparable to the 'raw' age predicted 
by WD cooling. 
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model-dependent estimate for the age of the disc, To, given 
by 



To = 7© + inuc + 5 



(2) 



In the chemical evolution models developed in CPT and 
CT (respectively without and with radial flows), there is no 
reference to the process of formation of the halo and there- 
fore no prediction of the time delay between the formation 
of the halo and the disc. In other words, our chemical evo- 
lution model does not allow a direct estimation of the age 
of the halo, which we denote Ta, nor of course the age of 
the universe, which we denote Tu, since our model places 
no constraints on the 'halo gestation period' for our galaxy 
- i.e. the time elapsed from the Big Bang to the formation 
of the halo. As we already remarked in section 1, however, 
for consistency we clearly require that the following chain 
of inequalities always holds: 



Tu >T G >T D > T s 



(3) 



In this paper we will derive constraints on the infall 
parameter, /3, of our chemical evolution model by compar- 
ing our model predictions for the age of the disc with other 
estimates of the age of the solar neighbourhood, the galac- 
tic halo and the universe, subject to the condition that the 
above inequalities should hold - in other words that our 
chemical evolution model should be consistent with other 
astrophysical and cosmological age determinations. 

In order to determine the abundance ratios of the ra- 
dioactive isotopes, we follow the same procedure as in CPT, 
with the simple difference that when SF starts, we solve the 
following equation for the evolution of the fraction of mass, 
Za, of a radioactive isotope A (Tinsley 1980) 



oIZa 
dt 



Pas - clZa 



Aa Za 



(4) 



where Aa is the rate of disintegration of the isotope, A, pa 
is its yield - assumed to be constant - E 3 is the gas surface 
density, s is the SFR and a — aof is the infall rate. We have 
used the model infall function 

21 

/ = 5/3exp(-fcr)[sech(/3t)(l + ltanh(/3t))]/(l + — ) (5) 

where the denominator is introduced as a normalisation fac- 
tor (to ensure that the final mass of the galaxy is not de- 
pendent on either /3 or I) and /3 is the infall parameter, as 
defined in CPT. We have chosen the value of the parameter 
oo such that the final mass of our model galaxy is in agree- 
ment with the observed range of values for our Galaxy. In 
conjunction with equation ^, we solve the equation for the 
gas surface density, 

dE 9 

—7 s - = a — s 
dt 



(6) 



As in CPT, we require that SF starts only when the 
Toomre stability criterion is marginally satisfied, i.e.:- 



Q 



7rGE 



(7) 



where k is the epicyclic frequency and c g the gas velocity 
dispersion. This threshhold on SF implies the delay, S, dis- 
cussed above - which is partly dependent on the thermal 
properties of the disc and the characteristics of the infall, 
amongst other physical properties (See CPT for a more de- 



tailed discussion). In the present paper we concentrate on 
the properties of the infall rate. 

We write equations ^ and ^ in their non-dimensional 
form 



dZA 

~dT 

and 

d:t 



y- f z a 



— Aa za 



dr 



= f-y 



(8) 



(9) 



where we have used the notation defined in CPT; in partic- 
ular here we have written Za = Paza ■ With this notation, 
the abundance ratio of an isotopic pair (A,B) with a pro- 
duction ratio Pa/Pb is 



N A = Pa za_ 
N b ~ Pb z b 



(10) 



Here we have used the definition of the yield as the 
ratio of newly synthesised elements returned instantaneously 
to the ISM, to the fraction of mass locked up in long-lived 
remnants: i.e. pa = Am u PA/(l — R), where R is the ejection 
rate and m u is the atomic mass unit. Furthermore we have 
denned A a = to A a as the dimensionless disintegration rate. 
Note that a fully consistent analysis should ideally consider a 
variable yield, especially at the early stages of the production 
of the radioactive elements (cf. Pagel 1993). 

In order to calculate the initial (i.e. at the time when 
SF begins) enrichment of a given isotope, A, we assume in 
our models that there existed a 'primordial' metallicity, zio, 
originating from earlier nucleosynthesis in the halo at an av- 
erage time of approximately t e Gyr before SF starts in the 
solar neighbourhood. This generic assumption can find some 
support in the observed high (« 0.5) [O/Fe] ratio for metal 
poor halo stars, indicating contamination from high mass 
stars formed during early SF in the galactic halo (Burkert 
1993). Under this assumption the initial metal enrichment 
would then be given by za% = z^o • exp( — AAt e ). To improve 
upon this approximation we would need to further improve 
the models of CPT by a detailed study of the early moments 
of galaxy formation, in order to be able to follow the chem- 
ical evolution of the halo and the halo/disc interaction. We 
will show in section 5, however, that the predicted abun- 
dance ratios at late times have negligible dependence on the 
adopted values of Zio and t e . 



2.2 Radial flows 

For the case where we have studied the effect of radial flows 
on predictions for the age of the disc, we have used the cal- 
culations developed previously in CT - with the only change 
in the metallicity equation, in which we have added the term 
— Aa Za E 9 to the right hand side to account for the decay 
of the radioactive element in the gas. As was shown in CT, 
after a certain time of evolution of the disc when SF is well 
underway, incorporating radial flows in the model will tend 
to produce a lower metallicity. Hence for a given (i.e. the 
observed) metallicity one would expect the disc to appear 
older under the influence of radial flows than in their ab- 
sence. This is because the fresh, metal-poor, gas flowing into 
the disc from the halo dilutes the heavy metals produced by 
SF. 
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The equation of evolution of radioactive elements is 



thus, 



v (dz A dz A \ 



clZa — Xa^oZa 



(11) 



Expressed in its non dimensional form this equation becomes 



2^ + Vr ^± 

dr dr 



y - f z A 



A.AZA 



(12) 



where f — r/ro, V r = v r /vor, vq t = ro/ta, with ro = 1 kpc 
and to = 1 Gyr. 

At this stage we need to treat carefully the solution of 
the above equation. When we include radial flows all regions 
of the disc are intercorrelated as star formation propagates 
outwards from the centre of the disc (see figure 1, where we 
plot the time, at which SF begins at radius, R*) and ide- 
ally we should account for this when modelling the initial 
metal enrichment of the solar neighbourhood. A fully consis- 
tent treatment to this problem is far from trivial, however, 
and - as we indicated previously - would require a detailed 
model for the dynamical evolution of the halo, both prior to 
the onset of SF in the disc and throughout the subsequent 
evolution of the disc. We can, nevertheless, place useful lim- 
its on the impact of including radial flows on our age predic- 
tions for the disc by considering the case where the initial 
metal enrichment is identically equal to zero. Since the effect 
of radial flows will always be to dilute the primordial metal- 
licity, Zio, setting the initial enrichment equal to zero when 
SF begins at each radius in the disc is simply the limiting 
case of this dilution process. 



3 RESULTS 

We now describe in detail the numerical results of solving 
the equations of our chemical evolution models outlined in 
section 3. Our solution involves the adoption of numerical 
values for the following three parameters: K\, Ki and Ks, 
which denote the fraction of energy put into the ISM by in- 
fall, SF and cloud-cloud collisions respectively. We have used 
the following numerical values (arguments for the choice of 
which were given in CPT): 



K x = 2 , K 2 = 25 , K 3 



10 



(13) 



The values of the (dimensionless) disintegration con- 
stants which we adopt are A A = 0.0495 for 232 Th, 0.985 for 
235 U and 0.1551 for 238 U (cf. equations | and|l|). Varying 
the remaining free parameters of our models, we found that 
the predicted abundance ratios are sensitive only to the in- 
fall parameter, f3. We thus chose to address quantitatively 
the sensitivity of our age predictions to /3 and assigned to 
all other model parameters the same numerical values as in 
'Model 1' of CPT (See Table 1 of CPT for further details). 

Figure 2 illustrates the dependence of SF on the infall 
parameter, plotting the SFR as a function of time for f3 = 
0.1, 0.4, 0.7 and 1.0, and it is clear that the value of /3 deter- 
mines not only the delay, S, before the onset of SF, but also 
the duration and amount of SF in the solar neighbourhood, 
at a given time. For /3 = 0.1 the onset of SF is delayed by 
S ~ 6 Gyr, and SF then proceeds as a slowly varying func- 
tion of time. For f3 = 1.0, on the other hand, the delay, S, is 
less than 1 Gyr, and is followed by a rapid burst of SF which 



has decayed to a few percent of its peak rate after only ~ 7 
Gyr. 

In order to address quantitatively this sensitivity to /3, 
we consider below a statistical analysis designed to obtain 
maximum likelihood estimates of f3, and thus constrain the 
law of SF, which 'best fits' our chemical evolution model 
to the observational age constraints in a manner consistent 
with the inequalities of equation ^. 

Among other tests to which one can subject our mod- 
els are their ability to solve the G-dwarf problem, to fit the 
age-metallicity relation and to predict the observed gas frac- 
tion in the solar neighbourhood. The relative success of the 
models in addressing these problems has already been dis- 
cussed in CPT. In figure 3 we illustrate the explicit depen- 
dence of the gas fraction, /i, on the adopted value of ft - 
associated with the SFR as shown in figure 2. We plot the 
model-predicted values of [t, as a function of the present age 
of the galactic disc, and compare them to the observed value 
of /iobs = 0.28 ± 0.09, taken from Gilmore, Kuijken & Wyse 
(1989), and denoted by the horizontal lines on figure 3. We 
also illustrate, shown on the upper ordinate of figure 3, the 
relationship between the model-predicted values of \i and the 
age of the universe - a lower limit for which is estimated by 
adding 1.5 Gyr to account for the time interval between the 
Big Bang and the formation of the galactic disc. 

We can see from the figure that the predicted value of {i 
matches the observed value at a wide range of different disc 
ages, depending on the value of /3, ranging from ~ 20 Gyr 
for j3 = 0.1 to ~ 3 Gyr for /3 = 1.0. N ote that for (3 > 0.4 the 
best fit of the model predictions to the observations is found 
for a low disc age of t < 6.5 Gyr - in good agreement with the 
age of the universe predicted for the standard inflation model 
with a fairly high (Ho ~ 70) Hubble constant, even if our 
adopted disc 'gestation period' of 1.5 Gyr were something 
of an underestimate. The shape of the curves and the width 
of the observational error band on jj, b B , however, mean that 
considerably higher disc ages are not strongly excluded - 
even for /3 = 1. 

We have computed the abundance ratios of the chrono- 
metric pairs 232 Th/ 238 U and 235 u/ 238 U with the respective 
production ratios 1.6±0.1 and 1.24±0.1 from Cowan, Thiele- 
mann & Truran (1991). Our computations show that the 
quoted uncertainty of ±0.1 on the production ratios intro- 
duce a spread of about ±1 Gyr on the predicted age, in the 
sense that the higher the production ratio the younger the 
disc and vice-versa. 

Figures 4a and 4b show the time evolution of the 
model predicted values of the abundance ratios of our two 
chronometers - again associated with the same SFR of fig- 
ure 2. As in figure 3, the observed abundance ratio and its 
la error band are denoted by the horizontal lines on each 
figure. These graphs have been obtained adopting a delay 
of t e — 3 Gyr between the period of nucleosynthesis in the 
halo - which we assume resulted in a 'primordial' metallic- 
ity of Zio = 0.25; see section 2 - and the onset of SF in the 
solar neighbourhood. Reasons for the choice of Zio = 0.25 
are discussed in CPT; we show in section 5, however, that 
the adopted values of both Zio and t e have negligible impact 
on our results. 

Note that, in a similar manner to figure 3, we estimate 
the time evolution of the production ratios in our model as 
a function of the age of the universe (as shown in the upper 



© 0000 RAS, MNRAS 000, 000-000 



Cosmochronology and the age of the galaxy 5 



ordinates of figures 4 and 5) by adding 6 Gyr to the age of 
the disc: a disc gestation period of 1.5 Gyr and a further 
4.5 Gyr representing the interval from the formation of the 
solar system to the present day - during which time the 
abundance ratios have been 'frozen'. As in the case of figure 
3, this second age scale is designed to be largely illustrative 
- simply in order to allow an approximate lower bound for 
the age of the universe to be estimated from our models. 

A quantitative examination of figures 4a and 4b re- 
veals several interesting features. Firstly, note that the ini- 
tial value of the 232r ph/ 238 U ratio predicted in our model lies 
within the la error band of the observed value of this ratio, 
and only deviates significantly from the observed ratio sev- 
eral Gyr after the onset of SF (the exact time depending on 
the value of (S) before reaching a minimum value and then 
increasing monotonically at later times. The 232 Th/ 238 U 
abundance ratio, on the other hand, is predicted in the 
model to have an initial value several standard deviations 
below the observed ratio, to increase rapidly after the onset 
of SF, to reach a maximum several Gyr later and then to 
decay monotonically at later times. 

A naive interpretation of these figures could then lead 
to the conclusion that our model predictions give a good 
fit to the observed abundance ratios for very young ages. 
For example, in the case of f3 = 0.4 we see from figure 4b 
that the ratio 232 Th/ 238 U rises sharply to intercept the ob- 
served abundance ratio when the age of the disc is only 
2 Gyr. At the corresponding age on figure 4a we see that 
the model predicted ratio of 232 Th/ 238 U has decreased only 
very slightly from its initial value, and still lies comfortably 
within the la error band. Thus, when /3 = 0.4, it is appar- 
ent that the observed and predicted abundance ratios are 
in good agreement for a disc which is only 2 Gyr old. The 
obvious flaw in this interpretation, however, is the fact that 
SF has begun only 1 Gyr before the formation of the solar 
system with these model parameters. In other words this 
apparently good fit to the observed abundance ratios does 
not allow a nearly sufficient duration for nucleosynthesis. 
We know, from the ages of WD in the solar neighbourhood, 
that significant amounts of SF had occured in our region 
of the galactic disc around 9 Gyr agoQ. Allowing 4.5 Gyr 
for the age of the solar system, this means that we should 
allow a minimum of ~ 4.5 Gyr to account for the duration 
of nucleosynthesis - i.e. the interval, t nuc , between the onset 
of SF in the solar neighbourhood and the formation of the 
solar system. In particular, therefore, any assessment of the 
'goodness of fit' between our model predictions and the ob- 
served abundance ratios must ensure that a sufficiently long 
period of nucleosynthesis is allowed. Similarly, a statistical 
analysis designed to estimate 'best fit' values of the infall 
parameter, (3, should exclude all cases in which the duration 
of nucleosynthesis is unreasonably short. We now describe 
such a statistical analysis. 

We begin by introducing some notation. Let U b s and 
Th bs denote the observed abundance ratios of 235 U / 23 8tj 
and 232 Th / 238 U respectively, and let Utruc and Th truc de- 
note the true values of these ratios. We define eu and e-rh, 
the observational errors on the measured values of the pro- 
duction ratios, as: 



t See also section 4 for a more detailed discussion 



Uo 



Utr 



£Th = Th D 



Thtruc 



Finally we define U mo dei(i, P) and Uth(t, P) as the present 
day values of 235 U / 238 U and 232 Th / 238 U predicted in 
our model as a function of the age of the universe, t, and 
the value of the infall parameter, /3. We wish to construct a 
joint likelihood distribution, C(t,/3), for t and j3 under the 
null hypothesis that our model is correct. In other words, 
C{t,f3) measures the probability that one would obtain the 
observed values of the abundance ratios, given that their 
true values are equal to the values predicted in our model - 
for given t and f3. Thus, under the null hypothesis, we have: 



eu = U bs — U n 



£Th 



= Th 



Th n 



We assume that eu and e-rh are normally distributed 
with zero mean and dispersion equal to the standard error 
of the respective abundance ratio, as reported in the liter- 
ature. We also assume that eu and £Th are uncorrelated - 
which is certainly not the case for the values of Utrue and 
Thtrue, but is likely to be a reasonable approximation for 
their observational errors. 

The joint likelihood function, £(t,(3), is then given by: 



C(t, f3)dtd(3 = Aexp 



Uo 



an 



1 / Th obs 

2 V 



Th n 



S(t, P)dtd/3 



(14) 



Here ctu and orh denote the dispersion of eu and eTh re- 
spectively, A is a normalisation constant and S(t, (3) is a 
selection function which excludes values of t and /3 which 
are physically unreasonable due to, e.g., too short a dura- 
tion for nucleosynthesis. 

It now follows from equation (jl4]) that the likelihood 
distribution of t conditional upon (3 is given by 



C(t\/3)dt = 



C(t,f3)dt 
J C(t,f3)dt 



(15) 



In figure 5 we plot this distribution for the illustrative 
values of f3 — 0.1, 0.4, 0.7, 1.0, as before. We assume a mini- 
mum duration, t n uc, of nucleosynthesis of 4.5 Gyr. Of course 
the minimum acceptable age for the disc at the formation 
of the solar system is in practice higher than 4.5 Gyr, since 
our model also predicts a j3 dependent time delay, 8, before 
the onset of SF in the solar neighbourhood, as we have al- 
ready seen in figure 2, and this limiting age causes the slight 
skewness in the likelihood distribution for ft = 0.7 and 1.0 
in figure 5. We see that for ft = 0.1 we require an age for 
the disc at solar system formation of between ~ 13 Gyr and 
~ 24 Gyr to match the observed abundance ratios with a 
non- negligible likelihood, while for (3 = 1.0 the likelihood 
is non-negligible over the much narrower range of ~ 5.3 — 7 
Gyr. Note that, while the likelihood distributions for /3 = 0.1 
and 0.4 are effectively disjoint, there is a significant overlap 
between the likelihood distribution for /3 = 0.4 and 0.7 and 
for /3 = 0.7 and 1.0. This reflects the fact that the likelihood 
distributions become increasingly 'pushed up' against the 
lower limit of t nuc > 4.5 Gyr as (3 increases. Another con- 
sequence of this is the fact that as t decreases, the value of 
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(3 is less well constrained by the observed abundance ratios, 
since there is a wider range of values of (3 for which C(t\/3) 
is comparable in magnitude. 

We can illustrate more precisely the constraints on (3 
as a function of t by considering the conditional likelihood, 
£(/3|t), given by 



C(l3\t)d/3 = 



C(t,p)dp 



(16) 



Figure 6 shows likelihood curves for £(/3|tdisc) for sev- 
eral values of tdisc, the age of the disc at the formation of the 
solar system. It is clear that the likelihood distribution does 
indeed become slightly narrower as the age of the disc in- 
creases. Using equation [l4| it is straightforward to construct 
a curve of the maximum likelihood estimate, jSmx (tdisc), of 
(3 as a function of tdisc. We plot this curve in figure 8. We 
can see that if tdisc lies in the range e.g. 8-12 Gyr, this 
corresponds to a maximum likelihood estimate of /3 lying in 
the approximate range 0.25 - 0.45. 

Of course the form of the likelihood curves plotted in 
figures 5-7 depends in part upon our adopted value for the 
minimum duration of nucleosynthesis - so far taken to be 
4.5 Gyr. Increasing t nuc to, e.g., 6 Gyr causes a further in- 
crease in the skewness of the likelihood distributions, C(t\/3) 
for large values of (3. The likelihood distributions conditional 
upon age, C(J3\t), are largely unaffected by the increase in 
tnuc, however, although of course there is now no likelihood 
curve for a disc age of 6 Gyr as this value is excluded by 
the higher value of t nuc because one must also allow for the 
f3 dependent time delay before the onset of SF in the so- 
lar neighbourhood. Figure 8 shows the maximum likelihood 
estimate of f3 as a function of tdisc, but now with t nuc con- 
strained to be at least 6 Gyr. We see that there now exists 
no acceptable solution for (3ml (t) for a disc younger than 
~ 7 Gyr. In other words, it is not possible for our model 
to build the disc quickly enough and still allow at least 6 
Gyr between the onset of SF and the formation of the Solar 
System. Adopting our estimate of 1.5 for the disc gestation 
period, this translates into a lower limit on the age of the 
universe, Tu, of 13 Gyr. If Tu lies in the interval ~ 13— 13.5 
Gyr, the effect of increasing the lower limit on t nuc to 6 Gyr 
is to force a sharp increase in the value of $m x(i), in or- 
der that SF begins early enough to give a sufficiently long 
duration for nucleosynthesis. Above an age of ~ 13.5 Gyr, 
however, /3m x(t) is identical to the value obtained with the 
smaller lower limit of t nuc > 4.5 Gyr, as shown in figure 7 - a 
fact which is easily verified analytically from differentiation 
of equation 

Following the method described in Hendry & Simmons 
(1990) we can construct confidence curves for the infall pa- 
rameter, f3, and the age of the disc, tdisc, when the solar 
system formed. From these curves we can derive confidence 
intervals for f3 given tdisc, or tdisc given (3. Examples of such 
curves are shown in figures 9a and 9b, where the confidence 
intervals are derived at 67% and 90%, and with an assumed 
lower limit of t nuc > 4.5. By drawing a vertical line at some 
fixed value of tdisc the points at which the line intercepts the 
two curves define the limits of the appropriate confidence in- 
terval for (3. 

We can see from figures 9a and 9b that as tdisc increases, 
the confidence intervals for f3 become somewhat narrower 
in width and the limits of the intervals are monotonically 



decreasing - i.e. the estimated value of f3 is both smaller 
and more tightly constrained for larger values of tdisc- We 
can see, moreover, that for tdi sc < 8 Gyr, the confidence 
intervals are pushed sharply to high values of j3. This effect 
is also dependent on the adopted lower limit for t nuc , and can 
be seen to be somewhat more pronounced for the confidence 
curves with t nuc > 6 Gyr, as shown in figures 10a and 10b. 
Clearly a high value of f3 is required if SF is to begin in the 
solar neighbourhood quickly enough to allow a disc age of 
less than 8 Gyr when the solar system is formed and at the 
same time allow a sufficient duration for nucleosynthesis - 
particularly with t nuc > 6 Gyr. As was the case with figures 7 
and 8, on the other hand, the adopted lower limit on t nuc has 
no effect on the shape of the confidence curves for tdisc > 10 
Gyr. 

In order to constrain (3 more precisely from this analysis 
- and to determine if the age predictions of our model can be 
made consistent with the inequalities of equation |-we now 
consider what are the constraints upon the age of the disc 
from estimates of white dwarf cooling and globular clusters. 



4 CONSTRAINTS ON THE AGE OF THE 

GALACTIC DISC FROM OTHER METHODS 

4.1 White dwarfs 

On the basis of the theory of cooling of white dwarfs (WD), 
developed by Mestel (1952), Schatzman (1958) and Mestel 
& Ruderman (1967), which gives a relation between their 
age and luminosity, one can estimate the age of the galactic 
disc by measuring the luminosity of the coolest white dwarfs 
(Schmidt, 1959). This method requires the adoption of an 
equation of state for white dwarfs, knowledge of their chemi- 
cal composition (i.e. the relative abundance of C, O, He and 
H) and the fluid - crystal phase transition, since the cores of 
white dwarfs are supposed to crystallize at low luminosity. 
We do no more than summarise some salient points from 
recent applications of this method. For a comprehensive re- 
view see e.g. D'antona & Mazzitelli (1990). 

Hotter WD cool more rapidly, therefore the space- 
density of WD is expected to increase monotonically with 
decreasing WD luminosity. On the other hand, because of 
the finite age of the galactic disc there should be a paucity 
of WD at low luminosity. This is confirmed by observations 
which show an abrupt fall off in the number of WD stars 
below log(L/Lo) = —4.5. Therefore, because the WD are 
considered to be the oldest stars in the disc, the age of the 
galactic disc can be inferred by computing the time required 
for them to cool to this luminosity, with a correction to ac- 
count for their main sequence (MS) lifetime. 

Winget et al. (1987) have fitted the observational data 
of WD luminosities with a model assuming a pure car- 
bon core and a constant birthrate of WD over the age 
of the Galaxy. They deduced an age of the galactic disc 
Twd = 9.3 ± 2 Gyr. They then adopted an estimate for the 
age of the universe of Tu = 10.3 ±2.2 Gyr, which they claim 
accounts for the "time between the Big Bang and the first 
appearance of stars in the Galactic disk" . It seems to us that 
their addition of only 1 Gyr to arrive at an estimate of Tu is 
at best somewhat arbitrary and at worst a considerable un- 
derestmate, given that their Twd is inferred only from WD 
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in the solar neighbourhood and hence only takes account of 
the history of SF at our location in the disc. Their value 
of Twd = 9.3 ± 2 Gyr does, however, provide a reasonable 
estimate of the age which we denote by Ts in equation ^ - 
the age of the solar neighbourhood. 

Iben & Laughlin (1989) considered in their model more 
ingredients, such as the IMF, time-variable SFR and the 
explicit influence of the lifetimes of MS progenitors. They 
found that the shape of the luminosity function is inde- 
pendent of any variation in the SFR, supporting thus the 
assumption of Winget et al. They deduced an age for the 
galactic disc of Twd =9 — 10 Gyr, consistent with the 
Winget et al. result. 

As with the Winget et al. study, we regard this result 
as a reliable estimate of Ts. Both determinations should be 
regarded as no more than a lower limit for the present age 
of the galactic disc, Td, and the age of the galaxy, Tg, for 
the following reasons: 

(i) disc stars in the solar neighbourhood form certainly 
after the halo and bulge stars with a time delay of a few 
Gyr 

(ii) not all the stars born in the disc stay there, especially 
the first generation of stars which inherit high random ve- 
locities from the gas out of which they formed. These stars 
therefore undergo more scattering, contributing thus to the 
WD paucity observed 

(iii) any peculiar phenomenon liberating energy (nuclear, 
gravitational etc.) will tend to prolong the cooling time, and 
therefore increase the age of the galactic disc obtained by 
this method. 

(iv) there is no reason why the evolution of WD should 
be independent of the chemical evolution of the disc, a fact 
pointed out in Pitts & Tayler (1992). In fact Pitts & Tayler 
showed, from their chemical evolution model, that the study 
of WD may be an underestimate of the true age: they de- 
duced an age ~ 1 Gyr higher than that predicted by Winget 
et al. 



4.2 Globular Clusters 

Globular clusters (GC) are known to be the oldest ob- 
served stellar systems in our Galaxy - and to some extent in 
the Magellanic Clouds. All studies have generally converged 
upon an age in the range Tgc = 13 — 15 ± 3 Gyr (Van- 
denBerg 1990; Renzini 1993), although the age probability 
distribution is highly skewed and renders ages below 12 Gyr 
very unlikely (Tayler 1986, Renzini 1993). 

The age of GC can be determined by comparing their 
observed H-R diagrams with theoretical H-R diagrams cal- 
culated as a function of time. Assuming that all the ob- 
served stars in the cluster have been formed on a time-scale 
which is short compared to the cluster age, and that all stars 
are formed from matter of the same composition, a direct 
measure of the age of the cluster can be obtained if the lu- 
minosity of the MS stars of given composition is known at 
the turnoff point (TO) - the position on the HR diagram 
which nearly coincides with the exhaustion of H at the cen- 
tre of evolving stars - by using the time-luminosity relation 
for the MS-TO (VandenBerg 1988). The two methods com- 
monly used to obtain the TO luminosity involve one of the 
following: 



(i) the direct MS fitting of observed colour-magnitude di- 
agrams to theoretical models 

(ii) the calibration of the magnitude difference between 
H-B stars and the MS turnoff, at the colour of the turnoff. 

The accuracy of either method depends upon how re- 
alistic is the theoretical H-R diagram used (Chiosi, Bertelli 
& Bressan 1992), as well as the extent of the uncertainties 
relating to the observed chemical compositions, reddening 
and distance moduli (Renzini 1993). It is these systematic 
effects which chiefly contribute to the quoted systematic un- 
certainty of ±3 Gyr in Renzini (1993). 

The GC 47 Tucanae has been the most studied and 
thus has the most reliably determined age (Hesser, Harris 
& VandenBerg 1987). Assuming an initial helium content 
of Y = 0.24 for the GC stars, its age was found to be 
T47TUC — 13.5 ± 2 Gyr. Note that the determination of 
the age of the GC is very sensitive to the uncertainty in the 
He abundances, whose minimum values should be bound 
by the primordial Big Bang abundance (Tayler 1986). Van- 
denBerg argues that an uncertainty as small as AY = ±0.02 
has considerable ramifications for GC ages. There is no clear 
age-metallicity relation for GC, however, even if the choice 
of the metallicity seems to have a great influence on the de- 
termination of the age - as shown in the case of the GC 
M92. Adopting a composition of Y = 0.24, [Fe/H] = -2.03 
and [O/Fe] — 0.70 for M92 yielded an estimated age of 
T 92 ~ 14 Gyr. Adopting [Fe/H] = 0, with the other ratios 
unchanged, gave T92 ~ 17 Gyr. One is forced to conclude 
that the absolute value of the age depends critically on the 
assumed chemical composition of the member stars. Hence 
any improvement in GC age determination relies, in part, 
on improved determinations of their chemical composition. 
Another important source of uncertainty is the stellar evo- 
lutionary model adopted. As an example, Vandenberg has 
shown that the inclusion of Helium diffusion - which has 
been shown to be an important process in Pop II stars - in 
the canonical models of stellar evolution can reduce the age 
of the GC by about 15% - pushing back towards the age 
predicted by WD. 

Chaboyer (1995) has also attempted to study the in- 
fluence on the GC age determination of the uncertainties 
in the input physics used in stellar evolution models, and 
concluded that the spread in 'acceptable' GC age determina- 
tions - given the systematic errors - could be at most 11 - 21 
Gyr. The two major sources of error identified in this study 
originate from the difficulty of modelling convection in stel- 
lar evolution models (e.g. changing the mixing length from 
1.5 to 3 added an uncertainty in the age of up to 10%), and 
from the uncertainties in distance determinations to GC. 
Currently, distance moduli cannot be determined to an ac- 
curacy of better than ±0.2 magnitudes, which contributes 
an uncertainty of up to 22% to the inferred GC age (Ren- 
zini 1993). Improvements in the calibration of each of the 
distance indicators commonly used - RR Lyrae variables, 
Subdwarfs and WD - are expected to come from the HST 
and the Hipparcos satellite in the near future, however. 

In summary, if we assume that GC are of a comparable 
age to that of the galaxy itself, then GC ages essentially 
provide a direct measure of Tg in equation ^[ Following 
Renzini (1993) and Chaboyer (1995), we adopt the estimate 
Tg = 13 — 15 ± 3 Gyr for the purposes of comparison with 
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our chemical evolution models - with the important proviso 
that ages of less than 12 Gyr are highly unlikely. 



5 DISCUSSION 

We now turn to the main question of this paper: how do 
the predictions of our chemical evolution models fit into the 
picture of other predictions for the age of the disc, and what 
constraints do they allow us to place on the value of /3? 

We have seen in figures 7 and 8 that the maximum 
likelihood value of (3 decreases monotonically as the age of 
the galactic disc increases. Suppose that the age of the disc 
at the formation of the Solar System is taken to be 9 Gyr 
- which would imply an estimate for Tu of ~ 15 Gyr, as- 
suming a reasonably short gestation period of 1.5 Gyr for 
the disc. This value would also be in agreement with the 
age predictions of globular clusters discussed above. In this 
case the maximum likelihood estimator of (3, /3ml = 0.39 
(cf. figure 8), and a 67% confidence interval for (3 is found 
from figure 9a to be 0.34 < (3 < 0.46 (assuming t nuc > 4.5 
Gyr; the confidence interval is very slightly wider if we re- 
quire tnuc > 6.0 Gyr). We can see from figure 2 that, for 
this maximum likelihood value of /?, SF began in the so- 
lar neighbourhood about 2 Gyr after the formation of the 
disc - i.e. about 11.5 Gyr ago, which is also perfectly consis- 
tent with the constraints on Ts from WD cooling. We can 
see, moreover, from figure 3 that f3 = 0.4 gives a reasonable 
fit to the observed gas fraction, for Tu = 13.5 Gyr. Recall 
also from CPT that this value of f3 was successfully able 
to fit the observed age-metallicity relation and to solve the 
G-dwarf problem. Thus, a value of (3 ~ 0.4 would seem to 
be favoured by our models - predicting as it does a disc age 
which is consistent with the acceptable ranges for the other 
disc age predictions which we have considered. 

The biggest difficulty in achieving consistency with the 
inequalities of equation H when (3 = 0.4 would, therefore, 
seem to be finding a viable cosmological model which allows 
Tu > 15 Gyr. This is not an easy task, however, if Ho lies in 
the range of 60 - 70 km s _1 Mpc" 1 , as a number of recent 
reliable determinations have reported (cf. Freedman et al. 
1994; Schmidt, Kirschner & Eastman 1994; Riess, Press & 
Kirschner 1995). Such a high value of Ho would certainly 
rule out the standard inflationary universe, with Qo = 1 
and Ao = 0. Adopting the lower bound for the matter den- 
sity, Q m > 0.35, deduced in Liddle et al. (1995), we find 
that we can have Tu ~ 15 Gyr provided Ao 7^ 0, an issue 
previously examined in Tayler (1986) in which a similar con- 
clusion was reached. Thus, it would appear that constraints 
on cosmological parameters still allow the consistency of all 
of the inequalities in equation ^ for a value of (3 ~ 0.4, al- 
though certainly the range of compatible values of Ho, fio 
and Ao is severely limited. In Hendry & Chamcham (1995, 
in preparation) we investigate in detail the current observa- 
tional status of these cosmological parameters, and address 
specifically the issue of how readily one can match the con- 
sistent age predictions for the galactic disc which we have 
discussed in this paper with the age of the universe predicted 
in cosmological models. 

We now consider the sensitivity of our results to the 
adopted value of the primordial metallicity, z^. Figure 11 
shows the 67% confidence curves for (3 and the age of the 



disc, assuming t nuc > 4.5 Gyr and now with z^ — 0. We can 
see immediately, on comparison with figure 9a, that for all 
values of j3 the change in the confidence curves is negliglible. 
Specifically, the confidence limits for the age of the disc at a 
given value of f3 are increased by ~ 0.1 Gyr (or equivalently 
the confidence limits for f3 at a given disc age are increased 
by ~ 0.01) when z^ is decreased from 0.25 to the limiting 
value of zero. Although the magnitude of the difference in 
the confidence limits between figures 9a and 11 is clearly 
negligible, it is worth noting that the sign of the difference 
is certainly in accordance with our physical understanding: 
when the primordial metallicity is reduced, one would expect 
to need a slightly longer time to build the observed heavy 
element abundance ratios, for a given value of j3. 

Note that figure 11 also demonstrates that our results 
are completely insensitive to the adopted value of t e '■ setting 
z^ — is numerically equivalent to letting t e tend to infinity. 
Although very large values of t e would, of course, not make 
physical sense - as in any case t e is constrained to be shorter 
than the time interval between the Big Bang and the onset 
of SF in the disc - it is, nevertheless, clear that the chosen 
value of t e is essentially unimportant to our analysis. 

Finally we consider the effect on our results of includ- 
ing radial flows in our chemical evolution models. As we 
explained in section 2, a completely self-consistent treat- 
ment of this problem would require a detailed model for the 
dynamical evolution of the halo. Instead, in this paper, we 
again consider the limiting case where Zio = everywhere 
in the disc at the onset of SF. Thus we explicitly include 
the effects of radial flows in the solution of equation |l^ only 
after SF begins, which is essentially equivalent to assuming 
that the inflow of fresh, metal-poor, gas from the halo prior 
to the onset of SF completely dilutes the metal content of 
the disc during that time. Figure 12 shows the 67% confi- 
dence curves for (3 and the age of the disc in this case, and 
assuming t lluc > 4.5 Gyr. We can see from this figure that, 
for a disc age of 9 GYr at the formation of the Solar System, 
the 67% confidence interval for (3 is now 0.40 < (3 < 0.52 - 
i.e. the upper and lower confidence limits for j3 are increased 
by 0.06. 

This small positive shift in the confidence interval for 
(3 is again in accordance with our physical expectations: in 
the radial flow model, the continual dilution of the disc with 
metal-poor gas from the halo would imply that we need a 
slightly higher value of (3 in order to build the observed heavy 
element abundances in a given time. The shift is clearly not 
a large effect, however, particularly since a more realistic 
model for the dynamical evolution of the halo prior to SF 
in the disc would always result in less dilution of the disc 
metallicity by radial flows, thus resulting in an even smaller 
shift in the confidence curves. Hence, it seems that a value 
of (3 close to, or very slightly larger than, 0.4 is the most 
favoured in order to achieve a consistent picture with the 
age predictions of white dwarfs and globular clusters. 



6 CONCLUSIONS 

In order to place constraints on the most sensitive param- 
eter of our chemical evolution models, we have compared 
in this paper the abundance ratios of heavy actinide pairs 
predicted by the threshhold galactic disc models of CPT 
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and CT with their currently observed abundance ratios in 
the Solar System. We have assumed that these nuclei are 
produced in direct proportion to the rate of star formation 
and calculated their abundances as a function of time, al- 
lowing for their decay. The values of the abundance ratios 
235 U/ 23S U and ^^/^X] at the time of formation of the 
solar system are well known and the aim of our calculations 
has been to discover at what disc age the observed ratios are 
attained. We have studied a series of disc models with all 
parameters fixed except one: the parameter /3 which deter- 
mines the rate of infall of matter to the disc. We also carried 
out some calculations varying the initial metallicity of the 
disc, produced by earlier stars in the halo, but found that 
this factor had negligible influence on the ages predicted by 
our models. Thus it seems that - at least in the context of 
our models - the heavy element abundances in the Galac- 
tic disc are determined essentially by the stellar component 
of the disc itself, and are independent of the history of the 
halo. 

In view of other evidence on the age of the oldest nearby 
stars, we have required that star formation began in the So- 
lar neighbourhood at least 9 Gyr ago, and have incorporated 
this constraint into a statistical analysis designed to deter- 
mine our maximum likelihood estimates f3 as a function of 
the current age of the galactic disc, and vice versa. 

Our model predictions showed that the deduced age of 
the Universe is a monotonically decreasing function of /3. 
A high value of (3 implies both rapid disc formation and 
a strong initial burst of star formation. A value of (3 ~ 0.4 
predicts the current age of the galactic disc to be ~ 13.5 Gyr, 
which is consistent with current estimates of globular cluster 
ages, and is also in agreement with the models considered 
in CPT which satisfy other constraints on the properties of 
the local galactic disc - such as the G-dwarf distribution, 
the age metallicity relation and the fraction of gas content. 
Lower values of (3 require a very large disc age to satisfy the 
cosmochronological contraints. 

We have also studied the effect on the age predictions for 
the galactic disc of incorporating radial flows in our chemical 
evolution models, following the treatment presented in CT. 
We simplified the analysis by assuming that the initial en- 
richment was reduced to zero everywhere in the disc in our 
models, We found that, even in this limiting case, the 67% 
confidence limits for (3 at a given disc age were increased 
by only ~ 0.05; we argued that any realistic, self-consistent 
model for the dynamical evolution of the halo and the disc 
under radial flows would result in less dilution of the disc 
metallicity, and thus result in a smaller positive shift in the 
confidence intervals for (3. 

We have also considered the consistency of the age de- 
termined by cosmochronology with other constraints on the 
age of the galaxy. The estimated disc age for f3 = 0.4 is larger 
than that deduced from white dwarf cooling, but within the 
uncertainties of that age. Our prediction is in any event per- 
fectly consistent with the ages of white dwarfs in the Solar 
neighbourhood since, as pointed out by Pitts and Tayler 
(1992), a delay in the onset of star formation means that 
the age deduced locally from white dwarfs is somewhat less 
than the true disc age. 

This leaves as the only major problem a reconciliation 
of the disc age with the age of the Universe obtained from 
current estimates of cosmological parameters. Of course this 



problem is essentially no different to that of reconciling glob- 
ular cluster ages with cosmology. If the value of the Hubble 
constant is indeed greater than 60kms _1 Mpc _1 , as recent 
observations suggest, then to achieve consistency between 
the age of the universe and both our predicted age of the 
disc and the age of globular clusters seems to require a non 
zero cosmological constant. If globular cluster ages could be 
reduced to as little as 11 Gyr - which has been proposed 
by Chaboyer (1995) as a very robust lower limit, and which 
would again make viable a cosmological model with Ho = 60 
and Ao = - then a higher value of (3 might be more appro- 
priate for our chemical evolution models. A value of (3 = 1 
can give a local disc age comparable to that of white dwarfs 
- indicating that star formation began in the Solar neigh- 
bourhood within 1 Gyr of the formation of the disc - and 
an age of the Universe of 11 - 12 Gyr. Such a model would 
imply the presence of luminous discs at much higher redshift 
than in the case of (3 = 0.4 and this would therefore provide 
a suitable observational test to discriminate between these 
two cases. We will investigate high (3 models, together with 
a more detailed study of the current status of cosmological 
parameters, in a forthcoming paper. 
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Table 1. The range of age predictions reported in the literature for the four methods of age determination considered in this paper 
Constraint Method Age, To, (Gyr) References 



age of the solar neighbourhood 
age of the galactic disc 
age of the galaxy 
age of the universe 



cooling of white dwarfs 9.3 ± 2 

cosmochronology 10 - 30 

globular clusters 13 - 15 ± 3 

cosmological parameters 8-30 



Winget et al. (1987) 
Cowan, Thielemann & Truran (1991) 
Rcnzini (1993), Chaboyer (1995) 
various, cf. van den Bergh (1994, 1995) 
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